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Abstract 

We present a systematic approach to thermodynamical identities 
and illustrate the power of these methods by displaying Mathematica 
notebooks to deal with a large variety of such identities. In concrete 
examples these can involve rather tedious, not to say impossible, com- 
putations when done by hand. 

1 Introduction 

The subject of thermodynamics is notoriously difficult for mathematicians. 
V.I. Arnold [Ar] famously put it in a nutshell as follows: 

Every mathematician knows that it is impossible to understand 
any elementary course in thermodynamics. 

He continues by explaining that 

the reason is that [the] thermodynamics is based on a rather 
complicated mathematical theory, on [the] contact geometry. 

In [Co] we presented an axiomatic approach to Gibbsian thermodynamics 
which avoided these difficulties and developed a systematic approach to the 
classical identities of thermodynamics. In this short article we wish to present 
a theory-free version of the latter topic with a view to making it directly 
accessible to readers without the theoretical framework presented in [Co]. 

The four Maxwell relations in thermodynamics, which are equivalent to 

the Jacobian identity -— — ^ — - = 1, are our starting point. Geometrically, 

d{p, V) 

they mean that the the corresponding map from the {p, \^)-plane into the 
(T, 5')-plane is area-preserving. 



1 



2 The notation 



In order to develop a systematic approach we prefer to employ a mathemati- 
cally neutral notation rather than the standard notation of thermodynamics. 
We begin by recalling the latter for the reader's convenience. 

Our starting point is the situation where we are given the temperature T 
and the entropy S as functions of the pressure p and the volume V. 

We use the following dictionary to jump between the purely mathematical 
notation and the thermodynamical one: u corresponds to T, v to S, p to x 
and V to y. For example, the thermodynamical equations 

T^pV, S= In (pV^) 

7-1 

of the ideal gas (we are omitting constants) corresponds to 

u{x,y) = xy, v{x,y) = — ^—-(\nx + j\ny). 

7 - 1 

It is a consequence of the Maxwell relations that one can define four 
energy type functions and we will discuss this in more detail below. 

3 Thermodynamical identities — the basic 
machinery 

We suppose that u and v are given as functions / and g oi x and y. Thus 
u — f{x, y), V = g{x, y). When J = 1, simple manipulations with differential 
forms proved the basic identities: 
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where we have highhghted the expressions which correspond to the Maxwell 
relations. 

In thermodynamic notation these are 

dT = fidp + f2dV, dS = gidp + g2dV 
dp = g2dT-g2dS, dV = -gidT + fidV 

dT ^ ^dS+—dp, dV ^ —dS-—dp 

92 92 92 92 

dT = ^dS-—dV, dp = —dS-^dV 
91 91 91 91 

dS = ^dT+ — dV, dp = —dT-^dV 

Jl Jl Jl 92 

dS = ^dT-^dp, dV ^ ^dT-^dp 

72 /2 /2 /2 

where we are using the key: u-H-T,v-H-S,x-<r^p, y-<r^V introduced above. 
In order to isolate the underlying patterns, we now use a numerical code. 
Thus 

u ^ 3^ T 
V ^ S 
a; — > 1 •<— p 
y ^ 2^ V. 

Partial derivatives will be denoted by triples in brackets. (3,1,2), for ex- 
ample, denotes in the neutral notation, in the thermodynamical 
one. In general, k) denotes the derivative of the variable z, regarded as 
a function of the j-th and k-ih. variable, with respect to the j-th variable. 

By reading off from the above list, we can express each partial derivative 
of the form {i, j, k) in terms of fi, f2, 9i, 92- For example: 

(3,1,2) = (3,2,l) = /2, (4, 1,2) = 51, (4, 2,1) = 52; 



(1,3,4) = 52, (2,3,4) = -51, (1,4,3) = -/2, (2,4,3) = A; 

and so on. We can then express any derivative (a, b, c) in terms of ones 
of the form (li, 1,2) or (e,2,l). Thus the four derivatives with x and v as 
independent variables can be computed:: 
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(3,4,1) 
(2,4,1) 
(3,1,4) 

(2,1,4) 



(3,2,1) 
(4,2,1) 

(4,2,1) 

(4,2,1) 
(4,1,2) 



(4,2,1) 

Then, as is standard in thermodynamics, we can introduce four energy 
functions E^^, E^"^, E"^^, E^^ such that dE"^^ ^udv-xdy, dE^^ = u dv+y dx, 
dE^^ = —v du + y dx, dE'^^ = —v du — xdy (the superfixes correspond to the 
independent variables — thus for E^^ these are x and u i.e. 1 and 3). 

We will discuss these in more detail below where the rationale of our 
notation will be explained. In terms of the classical notation: 

dE = T dS — pdV (energy) 

dE = —SdT — pdV (free energy) 

dG = -SdT + Vdp (Gibbs' potential) 

dH ^ TdS + Vdp (enthalpy), 

i.e., E^^ = G, = ^14 = H and = ^ 

If we arrange the energy functions in lexicographic order i.e. as E^^, 
^14^ ^23^ ^24 denote them by 5, 6, 7 and 8 in this order, then we can 
incorporate them into our system. For it follows from the definitions and 
simple substitutions that 

dE^^ = (y - vfi)dx - f2vdy 
dE^^ = {uf, + y)dx + uf2dy 
dE'^^ = -vfidx +{x - vf2)dy 
dE^"^ = ufidx + {uf2 - x)dy 

and so 

(5,1,2) = y-gh, (5,2,1) = -gh 
(6,1,2) = y + fgi, (6,2,1) = fg^ 
(7,1,2) = -gh, (7,2,1) = -x-gh 
(8,1,2) = fgr, (8,2,1) = -x + fg^. 

One of the potentially irritating features of the thermodynamical identi- 
ties is that many are related by a simple swapping of the variables while this 
is accompanied by changes of sign which seem at first sight to be random. 
The simplest example is displayed by the four Maxwell relations. We can 
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systematise such computations by introducing the symbol [a, b; c, d\ for the 
Jacobi determinant of the mapping (c, d) i-> (a, b) i.e. 



[a, b; c, d] = {a, c, d){b, d, c) — (a, d, c)(6, c, d). 

The determinant then takes care of the sign. 

uiu V I 

For example [3, 4; 1, 2] is the Jacobian ^, ' , and so is 1 (in this case this 

d{x,y) 

denotes the number 1), [3,2:4,11 is S ' and so is = — — = — ^ ' ' ^ . 

^' L ' ' ' J d{v,x) (4,2,1) 

Note that there are 1,680 such Jacobians. However, lest the reader despair, 
we will display below simple rules which allow us to express them all in terms 
of our primitive quantities (/ and g together with their partials and, of course, 
X and y). Later, we append a Mathematica notebook which computes all such 
expressions at the press of a button. 



4 Higher derivatives 

Some of the thermodynamical identities involve higher derivatives and we 
indicate briefly how to incorporate these into our scheme. We use the self- 
explanatory notation ((a, 6, c), d, e) for second derivatives. Thus ((3, 1, 2), 2, 1) 

/ Q f dE 

is just fi2- Note that this notation allows for such derivatives as I 7— I 



dT \dp 



which is ((8, 1, 2), 3, 4). Once again, we can express all such derivatives (there 
are now 18,816 of them) in terms of x, y, /, g and their partials (now up to 
the second order) using the chain rule. For 

((a,6,c),z,j) = ((a,6,c),l,2)(l,i,j) + ((a,6,c),2,l)(2,i,j) 

and (a, b, c) resp. (l.z, j) and (2, i, j) can be dealt with using the above tables. 

We display below a Mathematica notebook which computes all expres- 
sions of the form ((a, b, c),d, e). 



V/ s 



5 Derived quantities and thermodynamical 
identities 

The reason why there is a plethora of thermodynamical identities is simple. A 
large number of significant (and also insignificant) quantities can be expressed 
or defined as simple algebraic combinations of a very few (our primitive 
quantities x, y, f, g and their partials). Hence there are bound to be many 
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relationships between them. Our strategy to verify (or falsify) an identity is 
to use the above methods to express both sides in terms of these quantities 

and check whether they agree. 

Of course, there are myriads of such quantities and identities and we can 
only bring a sample. Thus we have 

the heat capacity at constant volume and 



dT), 



the heat capacity at constant pressure. In our formalism, cy = /(4, 3, 2) and 
Cp = /(4, 3, 1) and so, from our tables, 

f9i .92 
h J2 

Hence for the important quantities 7 = — and Cj, — cy we have 7 = 

^ cv im 

and Cp-cv^ f-TT' 
hJ2 

We illustrate our method by verifying the simple identity: 

Using the tables above, we can easily compute both sides in terms of our 

/ 

primitive expressions and get — — in each case. 

/1/2 

6 Computing (a, 6, c) and ((a, 6, c), ti, e) 

We can summarise these results in the following formulae: 

(a,l,2)(b,2,l)-(a,2,l)(b,l,2) 
^"'"'"'''J (c,l,2)(d,2,l)-(c,2,l)(d,l,2) 

and so 

^"^ ^'''^ (6,l,2)(c,2,l)-(6,2,l)(c,l,2) 
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which allow us to systematically compute any of the derivatives of the form 
(a, b, c) in terms of our primitives x, y, f, g, fi, /2, gi and g2- 
For the second derivatives we substitute 



(a,l,2)(c,2,l 



)-(a,2,l)(c,l,2) 



(6,l,2)(c,2,l 



)-(6,2,l)(c,l,2) 



into the formula 



(0, d, e) 



(0,l,2)(rf,2,l)-(0,2,l)(d,l,2) 
(rf,l,2)(e,2,l)-(d,2,l)(e,l,2) 



to compute ((a, b, c), d, e) in terms of our primitive terms (this time with the 
first and second derivatives of / and g). The advantage of these formulae is, 
of course, that one can write a simple programme to compute them. (It is 
always tacitly assumed in the above formula that the appropriate conditions 
which allows a use of the inverse function theorem hold). 

The formulae developed here suffice to create Mathematica notebooks 
to calculate any of the above quanitities both for the general case and for 
specific gas models and we display these below. 



We emphasise that the numerical code for the various thermodynamical 
quantities is a mere construct to facilitate their computation (ideally with 
the aid of suitable software) and the final goal is to express them all in terms 
of the basic quantities (x, y, f, g and the partials of the latter). It is then a 
routine matter to translate these into the standard terminology of thermody- 
namics if so required. For the convenience of the reader we give a dictionary 
of the relationships between them: 



Of course, p is pressure, V volume, T temperature, S entropy and W, F 
and E are free enthalpy, enthalpy, free energy and energy respectively. 



7 A notational survival kit 



1 2 3 4 5 6 7 8 

p V T S <^ W F E 

X y u V E^^ E^* £;23 E^\ 
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8 Mathematica notebooks 



8.1 A notebook for (a,6, c), [a,b;c,d] and ((a, 6, c), e) 

We start with a notebook which computes (a, b, c), ([o, b; c, d\ and ((a, b, c), d, e) 
for any combination of a, b and c between 1 and 8. 

ClearEffl, ff2, ff3, ff4, ff5, ff, gg] ; 
ffl[a_, 1, 2] = 
If [a == 1, 1, 
If [a == 2, 0, 
If [a == 3, D[ff [x, y] , x] , 
If [a == 4, D[gg[x, y] , x] , 
If [a == 5, y - gg[x, y] D[ ff[x, y] , x] , 
If [a == 6, y + D[gg[x, y] , x] ff[x, y] , 
If [a == 7, -gg[x, y] D[ff[x, y] , x] , 
If [a == 8, D[gg[x, y], x] ff[x, y]]]]]]]]]; 

ff2[a_, 2, 1] = 
If [a == 1, 0, 
If [a == 2, 1, 
If [a == 3, D[ff [x, y] , y] , 
If [a == 4, D[gg[x, y] , y] , 
If [a == 5, -gg[x, y] D [ff[x, y] , y] , 
If [a == 6, D[gg[x, y], y] ff[x, y] , 
If [a == 7, -X - gg[x, y] D[ff[x, y] , y] , 
If [a == 8, -X + ff[x, y] D[ gg[x, y] , y]]]]]]]]]; 
ff4[a_, b_, c_, 

d_] = (ffl[a, 1, 2] ff2[b, 2, 1] - 

ff2[a, 2, 1] ffl[b, 1, 2])/(ffl[c, 1, 2] ff2[d, 2, 1] - 
ff2[c, 2, 1] ffl[d, 1, 2]); 
ff3[a_, b_, c_] = ff4[a, c, b, c] ; 
ff5[a_, b_, c_, d_, 

e_] = (D[ff3[a, b, c] , x] ff2[d, 2, 1] - 

D[ff3[a, b, c] , y] ffl[d, 1, 2])/(ffl[d, 1, 2] ff2[e, 2, 1] - 
ff2[d, 2, 1] ffl[d, 1, 2]); 



Here we have defined three functions ff3, ff4 and ff5 of 3, 4 and 5 ar- 
guments, the variables having whole number values between 1 and 8. The 
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inputs a,b,c , a,b,c,d or a,b,c,d,e yield the values {a,b,c), [a,b;c,d\ and 
((a, b, c), d, e) respectively. In order to show that it can cope with the most 
complicated gas models we have used the following one which combines the 
Van der Waals gas and the Feynman gas (see [Co] for details). 

8.2 A synthesis 

A number of models of real gases can be subsumed in the following one. 

7, the adiabatic exponent, is a function of one variable and is a primitive 

of . This model has the double advantage that it incorporates the van 

7-1 

der Waals gas and also one introduced by Feynman which allows for the fact 
that in real gases, 7 depends on temperature. If the constants a and b vanish, 
we get what we call the Feynman gas (this model appears in [Fe] — see [Co] 
for a more detailed description). We get the van der Waals gas when 7 is 
constant and this specialises to the ideal gas when a and b vanish. 

We can now compute any of the above expressions for this model, using 
the above notebook. 



8.3 Thermodynamics identities and Grobner bases 

We can proceed further along the path to automatisation by invoking the 
concept and methods of Grobner bases. Rather than give a systematic treat- 
ment, we again display a corresponding notebook. 

The input below consists of two list of symbols which use an obvious 
code — x412, for example, denotes the partial derivative (4, 1,2). The second 
list is of the quantities which occur in the formulae we are interested in. 
They are regarded as variables and our setting is the ring of polynomials 
in these variables. This hst can be chosen at the discretion of the user 
(there are too many candidates to include all of them in one notebook). 
The first list is of the relationships between these symbols (each expression 
is to be thought of as being set equal to zero) beginning with the Maxwell 
relation followed by the formulae deduced above. In order to keep things 
simple, we have only used the relationships which we computed above. They 
can, of course, be enriched by triples which were computed by the above 
Mathematica programme. 

Clear [p, q] ; p = {x312 x421 - x321 x412 - 1, xl34 - x421, xl43 + x321. 
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x234 + x412, x243 - x312, x421 x341 - x321, x314 x421 - 1, 
x241 x421 - 1, x214 x421 + x412, x342 x412 - x312, x324 x412 + 1, 
x412 xl42 - 1, xl24 x412 + x421, x432 x312 - x412, x423 x312 - 1, 
xl32 x312 - 1, xl23 x312 + x321, x431 x321 - x421, x413 x321 + 1, 
x231 x321 - 1, x213 x321 + x312, x612 - x2 - x3 x412, 
x621 - x3 x421, x512 - x2 + x4 x312, x521 + x4 x321, 
x712 + x4 x312, x721 + xl + x4 x321, x812 - x3 x412, 
x821 + xl - x3 x421}; 
q = {xl, x2, x3, x4, xll2, xl21, x212, x221, x312, x321, x412, x421, 

xl34, xl43, x234, x243, x341, x314, x241, x214, x342, x324, xl42, 

xl24, x432, x423, xl32, xl23, x431, x413, x231, xl23, x512, x521, 

x612, x621, x712, x721, x812, x821}; 
GroebnerBasis [p, q] 

The output below is a (Grobner) basis for the ideal generated by the above 
list of relationships in our polynomial ring. Hence every entry is, when set 
equal to zero, a (new) identity. 

{x213 x621 + x712 + x213 x721 - x213 x821, x521 - x621 - x721 + x821, 
x512 - x612 - x712 + x812, 

x231 + x413, -x231"2 x621 + x431 x712 + x213 x431 x721 - x431 x812 - 

x213 x431 x821, -1 + xl23 x213, x621 + xl23 x712 + x721 - x821, 
xl23 x231~2 x621 + x431 x621 + xl23 x431 x812, xl32 + xl23 x231, 
xl23 x231 + x423, -xl23 x231'-2 - x431 + x432, x621 + xl24 x812, 
xl24 x231"2 - x431 + xl24 x213 x431, 

xl23 xl24 x231"2 - xl23 x431 + xl24 x431, 1 - xl24 x213 + xl42 x231, 
xl24 x231 + xl42 x431, xl42 + x324, 

x342 x712 + x213 x342 x721 + xl42'-2 x213 x812 - x342 x812 - 

x213 x342 x821, xl42 x213 + x231 x342, -xl24 x213 + x342 x431, 
x342 x621 - xl42"2 x812 + xl23 x342 x812, 
xl42"2 - xl23 x342 + xl24 x342, 

x214 x712 + x213 x214 x721 - x213 x812 - x213 x214 x821, 
x214 x621 + x812, -x231''2 - x213 x431 + x214 x431, -1 + 

xl24 x214, -xl42'-2 x213 x214 - x213 x342 + x214 x342, 
xl42 x214 + x241, xl42 x214 + x314, -xl42"2 x214 + x341 - x342, 
x243 x712 + x213 x243 x721 + xl42 x213 x812 - x243 x812 - 

x213 x243 x821, x213 + x231 x243, 
x243 x621 - xl42 x812 + xl23 x243 x812, xl42 - xl23 x243 + xl24 x243, 

xl42 x243 - x342, -xl42 x213 x214 - x213 x243 + x214 x243, -x231 + 

x234 + x243 x431, xl43 - xl23 x243, xl34 + xl23 x243 x431, 
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x421 + xl23 x243 x431, x231 + x412 - x243 x431, 
xl23 x243 + x321, -x243 + x312, 

x4 + x231 x621 + x231 x721 - x231 x821, x3 - xl42 x812, 
x2 - x612 + x812, xl - x621 + x821} 

The literature: The method of Grobner bases which is used in Mathe- 
matica is based on the so-called Buchberger algorithm (see [Bu]) and there 
is an extensive literature about it. 

We would also like to mention two works ([Br] and [Ja]) which are dedi- 
cated to the development of systematic approaches to thermodynamical iden- 
tities and which were useful in our approach. 
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